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ABSTRACT 

The origin of Galactic cosmic rays remains unconfirmed, but promising candidates 
for their sources are found in star-forming regions. We report a series of X-ray obser¬ 
vations, with Suzaku, toward the nearby star-forming region of Cygnus X. They aim at 
comparing diffuse X-ray emissions on and off the 7 -ray cocoon of hard cosmic rays re¬ 
vealed by Fermi LAT. After excluding point sources and small-scale structures and sub¬ 
tracting the non-X-ray and cosmic X-ray backgrounds, the 2-10 keV X-ray intensity dis¬ 
tribution is found to monotonically decrease with increasing Galactic latitude. This indi¬ 
cates that most of the extended emission detected by Suzaku originates from the Galac¬ 
tic ridge. In two observations, we derive upper limits of 3.4 x 10 -8 erg s _1 cm -2 sr ~ 4 
and 1.3 x 10 ~ 8 erg s _1 cm -2 sr ” 1 to X-ray emission in the 2-10 keV range from the 
7 -ray cocoon. These limits exclude the presence of cosmic-ray electrons with energies 
above about 50 TeV at a flux level capable of explaining the 7 -ray spectrum. They are 
consistent with the emission cut-off observed near a TeV in 7 rays. The properties of 
Galactic-ridge and local diffuse X-rays are also discussed. 

Subject headings: ISM: cosmic-rays — X-rays: diffuse background — X-rays: ISM - 
gamma rays: general 


1. Introduction 

An important question of modern astrophysics is the origin and propagation of Galactic cosmic 
rays (GCRs), which are charged particles with relativistic energies (up to 10 15 ~ 16 eV) diffusing in 
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and around the Milky Way. At low energy, they directly affect the chemistry and thermodynamics 
of interstellar matter through ionization, heating, and pressure. Supernova remnants (SNRs) are 
widely considered as the most plausible source s of GCRs, because they are e nergetic and numerous 
enough to maintain the power of GCRs (e.g., Ginzburg & Svrovatskiil 1964J). Because the massive 
OB stars (the progenitors of core-collapse supernovae) are born in clusters and live shortly, su¬ 
pernova expl osions tend to c l uster in space (within a few parsecs) and in time (within a few 10 5 
years) (IHierdon V Lingenfelterl l2005l. Therefore, GCRs are expected to be accelerated and injected 
into the interstellar space, not only by the ensembles of individual SNRs, but also b y overlapp ing 
shock s from SNRs and massive stellar wind s (called superbubbles. iTen orio-Tagle &; B odenheimer 


1988 1 created around OB associations (e.g., Bvkov &; Fleishman 19921 : 


Parizot et al 


2004). How¬ 


ever, this scenario about GCR origin requires observational confirmation. 


Recently, the Large Area Telescope (LAT) on the Fermi Gamma-Ray Space Telescope ([Atwood et al 


200 91) has revealed an extended source of hard, multi-GeV 7 rays, called the ’’ Cygnus cocoon 


(Ackermann et al. 

201 

l)i in 

19521: 

Uvaniker et al. 

1001 ) 


1.4 kpc flRygl et al. 2Q1_2|) • The extended emission 


is detected with high confidence above the interstellar 7 -ray background and after subtraction of 
known point sources, at energies above 1 GeV. The emission morphology corresponds to the region 
bounded by the ionization fronts powered by the numerous OB stars present in the region. The mor¬ 
phology and the absence of marked spectral variations across the cocoon imply an interstellar origin 
rather than a superposition of unresolved 7 -ray sources. The observed 7 -ray spectrum cannot be 
explained by the GCR proton or electron spectra measured at the Earth and in the local interstellar 
medium. The latter is often referred to as the Local Interstellar Spectrum (LIS). An amplification 
factor of (1.5-2) x (E p /10 GeV ) 0,3 of the LIS proton spectrum, or of ~ 60 x (E e /10 GeV ) 0 ' 5 of the 
LIS electron spectrum, is required to explain the cocoon emission (where E p and E e respectively 
denote the energy of the CR proton and electron in the cocoon). The intense and hard 7 -ray 
spectrum thus indicates the presence of freshly-accelerated CRs in the cocoon. 

Fermi -LAT data alone, however, cannot constrain the dominant type of radiating particles 
(protons or electrons), nor the maximum energy of the electrons. For this purpose, 7 -ray observa¬ 
tions are usefully complemented by X-ray observations, which are sensitive to CR electrons at TeV 
energies. We have therefore conduc ted deep X-ray obser vations of the 7 -ray cocoon r egion using the 


X-ray Imaging Spectrometer, XIS (jKovama et al.l 120071 ) on board Suzaku satellite (jMitsuda et al 


2007 ). which is very sensitive to extended X-ray emission. If the cocoon 7 -ray emission is produced 
by interactions of CR electrons with the interstellar radiation field and the particles accelerated 
beyond 10 TeV in energy have not yet cooled or escaped the region, we expect to find extended 
non-thermal X-ray em ission due to synchrotr on radiation in the large magnetic field of ~ 20 //G 
inferred in the region (Acker mann et ah 120111. 
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Observations and Data Reduction 


To constrain the CR properties in the cocoon, we have performed a series of X-ray observations 
of the Cygnus region. Two “source” observations were pointed at positions with strong GeV 
cocoon emission, but free of known, bright X-ray sources, such as the SNR 7 Cygni and the binary 
Cygnus X-3. Because a potential X-ray counterpart to the cocoon is likely to fill the whole field 
of view of the XIS instrument, we have also conducted two “background” observations to estimate 
the backgrounds i n the Cygnus direc t ion, in particular the intensity of the Gal actic ridge X-ray 
emission (GRXE; Worrall et al. 1982 : Warwick et al. 1985 : Kovama et al. 1986 ). The GRXE is 
an apparently diffuse feature along the Galactic plane and its origin would be a mixture of truly 


diffuse hot plasma and num erous dim point sources (e.g., lEbisawa et al .1 120051 : iRevnivtsev et al 


20061 : Uchiyama et al. 2013|). The positions of the four observations are overlaid on the 7 -ray count 
map of the cocoon in Figure 1. 

T he observations were ca rried out using the XIS on the focal plane of the X-Ray Telescope 
(XRT: iSerlemitsos et al.l 2007) on Suzaku. The XIS consists of two front-illuminated (FI) X-ray 
CCDs (XISO and 3) 0 and one backside-illuminated (BI) X-ray CCD (XIS1). The combined XIS 
and XRT system is sensitive within the energy range 0.3-12 keV. Although its angular resolution is 
moderate (h alf-power diameter ~ 2'), the XIS+XR T system provides a low and stable instrumental 
background (IMitsuda et al.ll2007l : lTawa et al.ll2008l ): therefore it is suitable for the study of extended 
emission with low surface brightness. The XIS was operated in the normal clocking full-frame mode, 
and data were recorded in the 3 x 3 or 5 x 5 editing mode. Data were analyzed using the HEADAS 
6.15.1 software with the calibration database released on July 1, 2014. We reprocessed data using 
aepipeline to take the latest calibration into account. We analyzed the so-called cleaned events, 
which had passed the following standard event selection criteria: a) only ASCA-grade 0, 2, 3, 4, 
and 6 events were accumulated with hot and flickering pixels removed; b) more than 436 s had 
elapsed since passing through the South Atlantic Anomaly; and c) the pointing direction was at 
least 5° and 20° above the rim of the Earth during night and day, respectively. To further reduce 
the non-X-ray background (NXB), we also required that d) the geomagnetic cutoff rigidity exceeded 
6 GV. Details concerning the observation and net exposures of the screened events are summarized 
in Table 1. 


1 Because of an anomaly occurred in 2006 November, the operation of another FI sensor, XIS2, has been terminated. 
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Fig. 1.— Positions of X-ray observations overlaid on the 7 -r ay residual coun t map of the cocoon 
above 10 GeV, smoothed with a a = 0.°25 Gaussian kernel (jAckermann e_t ah 2011;. The unit of 
the 7 -ray map is counts per pixel after the smoothing was applied. The coordinates are in Galactic 
longitude and latitude. Source and background observations are indicated by thick red and thin 
green boxes, respectively. See Table 1 for details. 


Table 1: Summary of observations 


Region 

Pointing* 1 

1 (deg) b (deg) 

Observation date 

Net exposure 
(ks) 

Source 1 

79.25 

1.49 

2012 Nov 18 

43.3 

Source 2 

78.99 

1.86 

2013 Dec 07 

42.0 

Background 1 

78.00 

0.74 

2012 Nov 19 

19.9 

Background 2 

80.50 

2.24 

2012 Nov 19 

25.6 


l and b are Galactic longitude and latitude of the center of the XIS field-of-view, respectively. 
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3. Data Analysis and Results 
3.1. X-ray Images 


We have extracted X-ray images from the two FI CCDs (XISO and XIS3). These cameras have 
better imaging quality than the BI CCD (XIS1) due to their lower instrumental background. We 
have defined the soft and hard bands as 0.7-2 keV and 2-10 keV, respectively, and have excluded the 
corners of the CCD chips illuminated by the 55 Fe calibration sources. We have then estimated the 
NXB cont ri bution from the night Earth data and subtracted it from the images using xisnxbgen 
(ITawa et al. 2008). Vignetting was then corrected by dividing the soft- and hard-band images 
by flat sky i mages sim u lated at 1 keV and 4 keV, respectively, by using the XRT+XIS simulator 
xissim (jls hisaki et al. 120071 1. In the flat image simulations, we have assumed a uniform intensity 
of 1 photon s -1 cm -2 sr -1 , thus the approximate unit of the obtained vignetting-corrected image 


is photons s _i cm -2 sr . We have combined two FI CCD images and applied smoothing with 
a Gaussian kernel of a = 0128 for visualization. Several point sources and small-scale structures 
are recognizable in the final images (see Figure 2). Moreover, we also observe strong, apparently 
extended emission in the soft-band images of Source 2 and Background 2 observations. In the 
hard-band images, the apparently extended emission is more intense in Source 1 and Background 1 
than in the other two regions. 


3.2. Extended Emission Spectra 


We have extracted extended emission based on the vignetting-corrected images (Figure 2). 
We have first identified point sources and small-scale structures in each of the soft- and hard-band 
images by eye and have delineated them by 1'5 radius circular regions (green circles in Figure 2). 
These structures have been designated as srcl, src2, etc. in each field, and events in those circles 
have been removecd. Spectral analysis of these sources and structures is described in § 3.2 and 
results are summarized in Table 3. An 80 pixel x 1024 pixel rectangular region and an area with 
(occasionally) significantly high background at the corner of the chip are unusable in XlSlIE For 
simplicity, these regions have been excluded from the sky projections for all XIS camera images. 
We have then accumulated the extended X-ray emission spectra by integrating the events over the 
whole XIS field of view (~ 17f8 x 17. / 8), and analyzed the spectra for diffuse X-ray emission with 
ancillary response files (ARFs) calculated by xissimarf gen ( Ishisaki et al12007 1. assuming uniform 
brightness over a 20' radius. The loss of effective area due to the exclusion of point sources and 


2 Although the choice of sources and structures is not unique, the effect on the estimated extended emission is small. 

For example, even if we included src3 (a moderately bright source located close to the center of the field of view) of 
Source 1 observation, the 2-10 keV intensity of the extended emission was increased by ~ 0.3 x 10~ 8 erg s _1 cm -2 sr^ 1 
which is less than half of the uncertainty examined in § 4.1 and shown in Figure 5 

3 See ’’The Suzaku Data Reduction Guide” of v5.0 (http: //www .heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc ) 
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unusable area described above has been taken into account in calculating the ARFs. The response 
matrix files (RMFs) has been calculated with xisrmfgen, and the NXB spectrum int egrated over 


the s ame region as that for the extended emission has been estimated by xisnxbgen (Tawa et al 


2008) and subtracted from the source spectrum. In our observations, the energy resolu tion degraded 


by particle radiation was restored by a spaced-row charge-injection technique (e.g., lUchivama et al 


2009al ) with a 6 keV equivalent charge applied on orbit. The amount of the charge injected affects 
the NXB level of the XIS, and was taken into account in calculationg the NXB. 

The NXB-subtracted X-ray spectrum is expected to be a sum of the cosmic X-ray background 
(CXB), t he GR XE (plus local diffuse X-ray emission), and the possible X-rays from the cocoon. 
Following I Kush ino et al. (2003), the CXB has been modelled as a fixed power law with a photon in¬ 
dex of 1.4 and a normalization of 9.42 x 10” 1 photons s” 1 end’ keV” 1 at 1 keV, or a 2-10 keV flux 


of 6.20 x 10 12 erg s 1 cm 2 , corresponding to an intensity of 5.85 x 10 


erg s 1 cm 2 sr 1 inte¬ 


grated throughout the 20' radius. Adopting the optically-thin approximation, from [K alber la et al 


(1200 51) we have estimated the neutral hydrogen column density toward the regions studied here 
1.3 x 10 22 cm -2 . However, the opacity of the Hi 21 cm line is non-negligible in the Galac- 


as ~ 


tic plane (e.g., 

Ackermann et al. 

2012g 

), and contribution from molecular clouds and dark neu- 

tral medium (e.g., 

Grenier et al. 

20od 

should be taken into account. We therefore estimated 

total neutral hydrogen column density from y-ray data in these regions based on the analysis by 


Ackermann et a l. (12012bi ). We adopted spin temperature of 250 K as a baseline model and ob¬ 
tained the total column density Y(H) of (2.5-3.1) x 10 22 cm” 2 , giving the 2-10 keV intensity of 
(4.8-5.0) x 10” 8 erg s” 1 cm” 2 sr” 1 (see Table 2). The uncertainty due to the assumed spin temper¬ 
ature is_evaluatedAn§Tl^_^ased_on_greyious_studies_ofthe_GRXEand local diffuse X-ray emission 


(e.g. JKaneda et al .11 199 71 : llJchivama et al.ll2009bl : lKataoka et al.ll2013l b we modeled the GRXE (plus 


local diffuse X-ray emission) as vapecl + wabs2 x apec2 + wabs 3 x apec3 in XSPEC, where apec (and 


vapec) and wabs are thin-thermal pla sma models by I Smith et al.1 (120011 ) and interstellar absorp 


tion by [Morrison & MCammo n (1983), respectively. We have allowed the relative normalizations 


among the three XIS sensors to vary. Among three components, the first term (low-temperature 
plasma model) represents the unabsorbed L ocal B ubble emission and/or contamination from the 


Solar-Wind Charge eXchange (SWCX) (e.g.. iFuiimoto et alJl2007l ). and we hav e fixed the plasma 


temperatur e a nd ab undance to kT = 0.1 keV and Z = Z®, respectively (see, e.g. JSmith et al.ll2007l : 


Yoshino et al. 20091). The second and the third terms (mid-temperature and high-temperature 


plasma models) correspond to t he classical soft-temperature and hard-temperature emissions of 


the GRXE, respectively (e.g., Kaneda et al. 19971)• In all four observations, we have found that 


our low-temperature plasma model is not able to reproduce the data in the energy range of 0.5— 
0.6 keV, because the observed energy of the emission line is lower than the O vn resonance line 
energy (0.574 keV). We have therefore set the Oxygen abundance to 0 and added a gaussian line 
at ~ 0.56 keV with the line center energy and normalization allowed to vary. The line is likely to 
be due to a charge-exchange process, as described in § 4.2. We have also added lines at ~ 0.65 keV 
(to reproduce the line of O vm) and at 6.40 keV (to reproduce Fe K a line) if the presence of the 
line was significant at the 99% confidence from an F-test. Our CXB+GRXE (plus local diffuse 
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X-ray emission) model successfully reproduces the Source 1 observation data (reduced chi-square 
X 2 /d.o.f. = 538.1/501). As shown in Figure 3, low-, middle- and high-temperature plasma models 
dominate the spectrum below ~ 0.7 keV, from ~ 0.7 to ~ 2 keV and above ~ 2 keV, respectively. 
Above 2 keV, the contributions from the CXB (fixed to a model described above) and the high- 
temperature plasma model are comparable. These two components are degenerate to a large extent 
in spectral shape, and the effect of the CXB uncertainty on the remaining hard-band intensity is 
described in § 4.1. The results from the fits in the four regions are summarized in Table 2 and 
Figure 4. 


3.3. Spectra of Individual Point Sources and Small-scale Structures 

To study point sources and small-scale structures, we have integrated the events within the 
If5 radius circular regions shown in Figure 2. If the If5 radius was too large and contaminated 
from nearby sources (e.g., srcl and src2 in Source 1 observation) or if the source was close to 
the edge of the CCD chip (e.g., src6 in Source 1 observation), the radius has been reduced to 1'. 
The RMF and ARF have been calculated by xisrmfgen and xisarfgen, respectively, assuming a 
point-source response in the latter. For the spectral fitting, we have used an absorbed power-law 
model (wabs x pow in XSPEC) or an absorbed thin-thermal plasma model (wabs x apec in XSPEC). 
We have first tried a power-law model and employed a thin-thermal plasma model if the power- 
law model was inadequate; i.e., if we observed large residuals in the fit or if the best-fit photon 
index were excessively large (T > 3). Src3 in Source 1 observation was not well fitted by a single 
component model, so we have fitted its spectrum using a two-component model (apec+pow with 
the common absorption). The fits results are summarized in Table 3. The majority of sources 
have a hard power-law spectrum with a photon index T of < 2.5 and relatively large absorption 
(A^(H) > 2 x 10 22 cm -2 ). These sources (src2, src4 and src5 in Source 1 observation, srcl and src2 in 
Source 2 observation, srcl in Background 1 observation and srcl in Back g rou nd 2 observation) are 
most likely background active-galactic nuclei. According to Tozzi et ajJ (2001) who compiled the 
log N - log S relation of X-ray sources, the expected number of active-galactic nuclei with 2-10 keV 
flux > 0.5 x 10 -13 erg s -1 cm -2 is ~ 30/deg 2 . This implies 2-3 uncatalogued active-galactic nuclei 
in each XIS field-of-view, which is close to the number of sources we found. Possible counterparts, 
identified from the SIMBAE0 and NEe| databases within V of the source position, are also listed 
in Table 3. 


J http://simbad.u-strasbg.fr/simbad/sim-fcoo 
J http://ned.ipac.caltech.edu/ 






Table 2: Best-fit spectral parameters of extended emission 



Sourcel 

Source2 

Backgroundl 

Background2 

Ei (keV) 

0.548 ± 0.004 

0.545 ± 0.006 

0.551 ±0.007 

0.554l 9 ;^ 

Normi 

90 o+3.7 
^•°-3.9 

24.4 ±3.1 

22 . 11 ^ 

14 -9l+o 

E 2 ( keV) 

0.653 (fixed) a 

0-648lffi 

0 673 +0 ’ 031 

U.O 1 >J_o 035 

- 

Norni 2 

9 q-)- 0.8 

z - J -0.9 

5.3 ± 1.1 

2 0 + °- 9 
z ’ u - 0.8 

- 

E 3 ( keV) 

6.40(fixed) 

- 

- 

- 

Nor m 3 

n 19 + 0.09 
U ’ 1Z -0.10 

- 

- 

- 

iV(H) low (10 22 cm" 2 ) 

0 (fixed) 

0(fixed) 

0(fixed) 

0 (fixed) 

kT low ( keV) 

0.1 (fixed) 

0.1 (fixed) 

O.l(fixed) 

0.1 (fixed) 

Mow {Zq) 

1 (fixed) 

1 (fixed) 

1 (fixed) 

1 (fixed) 

EM\ ow 

103 ± 18 


QQ+ 23 

^^-24 

139±|| 

N( H) mid (10 22 cm” 2 ) 

n on+0.07 
U - OJ -0.09 

0-46 

n 70 + 0.11 
u -' J -0.14 

0 (< 0 . 02 ) 

kT mid ( keV) 

fj yq+ 0.02 

u -'°-0.03 

0.62 ± 0.02 

0 62 +0,n 
u - dz -0.04 

0 60 + °' 02 
U ’ DU - 0.01 

Mnid(M-)) 

0-13iao4 

0 29 -1 " 0,27 
u - zy -o.io 

0 27" 1-1,25 
u - z '-0.15 

0-37i HI 

EM m i d 

256 

249+64 

zq:Z -98 

( 20 ±? 6 ) x 10 

73l3i 

^(H) high ( 10 22 cm” 2 ) 

0 1 + 1.4 

°' i -1.5 

1.9(< 3.0) b 

2 6 +0 ' 8 
z -°-i.o 

2 0+°' 7 

Z-U — 1.4 

feThigh(keV) 

0 1+2.4 
°' i -0.9 

2 - 6 +6 

2 1 + 0.8 

Z - 4 -0.4 

1 c+1.7 
i -°-0.3 

Miigh(-Z©) 

0 - 36 ^ai7 

0-28±^ 

n iq+°- 16 
u - iy -0.l4 

0.4(< 1.1) 

EM high 

106±|| 

7611 

217+ 78 
zi '-65 

IO 3 I 37 

/low (0.5-10 keV) c 

2.46 x 10 ” 8 

2.86 x 10” 8 

2.38 x 10 ” 8 

1.55 x 10 ” 8 

/mid(0.5-10 keV) c 

5.66 x 10 ” 8 

6.38 x 10 -8 

3.01 x 10 - 8 

8.42 x 10 ” 8 

/hi g h(0.5-10 keV) c 

4.88 x 10 ” 8 

3.15 x 10 ” 8 

7.20 x 10 ” 8 

2.39 x 10 ” 8 

1V(H) cxb (10 22 cm” 2 ) 

2.92(fixed) 

2.47(fixed) 

3.07(fixed) 

2.67(fixed) 

/cxb (0.5-2 keV) c 

0.12 x 10” 8 (fixed) 

0.16 x 10” 8 (fixed) 

0.11 x 10 ” 8 (fixed) 

0.14 x 10” 8 (fixed) 

/cxb(2-10 keV) c 

4.85 x 10 ” 8 (fixed) 

4.97 x 10 ” 8 (fixed) 

4.81 x 10” 8 (fixed) 

4.92 x 10” 8 (fixed) 

/o.5-2 C 

7.69 x 10 ” 8 

9.22 x 10 -8 

5.72 x 10 ” 8 

10.44 x 10 -8 

/ 2 -10 C 

10.29 x 10 ” 8 

8.33 x 10 -8 

11.79 x 10 -8 

6.98 x 10 ” 8 

X 2 /d.o.f. 

538.1/501 

565.7/521 

320.1/252 

370.5/336 


E and Norm are line center energy (keV) and the intensity (photons s 1 cm 2 sr 1 ) of gaussian, respec¬ 
tively. iV(H)i ow /fcTi OW Mi OW j E low; ;V(ff)mid /A.'/ in id/+i mid j d/mid and .V (TT j h i n 11 / /d/1 j n 11 / .4 hi g h / if iff high are absorp¬ 
tion/temperature/abundance/emission measure of the low-, middle- and high-temperature plasma models, respec¬ 
tively. The emission measures are given as the value integrated over the line-of-sight, ± f n e n^_ds (where n e and hh 
are the electron and hydrogen density, respectively) in the unit of 10 14 cm” 5 sr -1 . /i ow , /mid and /high are absorption- 
uncorrected intensity of low-, middle- and high-temperature components, respectively, where associated line fluxes are 
also included. 1V(H)cxb an< ^ /cxb are the assumed column density and absorption-uncorrected intensity of the CXB 
model, respectively. / 0 . 5-2 and / 2-10 are intensity in 0.5-2 and 2-10 keV including CXB contribution, respectively. 

“The line center energy was not well determined and thus fixed at the energy of O VIII K a line. 

6 The Lower limit was set at 1.0 x 10 22 cm -2 , otherwise the fit gives unphysical results. (Amid ~ 5, ^4high ~ 0 and 
N(H)hig h - 0) 

“The unit is given in erg s -1 cm -2 sr” 1 . 
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Fig. 2.— From top to bottom: X-ray intensity maps in Galactic coordinates of Source 1 observation 
(panels a), Source 2 observation (panels b), Background 1 observation (panels c), and Background 2 
observation (panels d). Left and right panels show the soft- (0.7-2 keV) and hard-band (2-10 keV) 
images, respectively. 
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Fig. 3.— Data and modeled spectra of extended X-ray emission of Source 1 region, taken by 
XIS1. The spectrum is fitted by three-temperature thin-thermal plasmas (GRXE plus local diffuse 
X-rays) and a fixed power law (CXB). Low-, middle-, and high-temperature models are indicated 
by red, green, and blue dotted lines, respectively, and the CXB contribution is shown by the dotted 
purple line. The bottom panel shows the residuals. 
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Energy (keV) Energy (keV) 

Fig. 4.— X-ray spectra of extended emissions integrated over the XIS field-of-view of Source 1 
observation (top left), Source 2 observation (top right), Background 1 observation (bottom left), 
and Background 2 observation (bottom right). The XISO, 1, and 3 data are indicated by black, 
red, and green crosses, respectively. Solid lines are the best-fit models. The bottom panels show 
the residuals. 




















































































Table 3: Spectra of point sources and small-scale structures 


Region/source 

Nh 

r 

kT 

A 

/o.5-2 

/ 2-10 

X 2 /d.o.f. 

possible 


10 22 cm ' 2 


keV 


erg s — 1 cm ^ 

erg s 1 cm ^ 


counterpart 3- 

Source 1 









srcl 

0.03(< 0.09) 

- 

0 94 + 0-08 
u - y ^-0.14 

0.07 ± 0.03 

1.81 X 10 -13 

0.29 X 10 “ 13 

71.6/64 

IRXS J202725.2+405428 

src 2 

5 1 + 5 - 2 
0 , 1 -2.7 

2 3 + 1-5 
z - 1.0 

- 

- 

0.03 X 10 -13 

1.52 x 10 “ 13 

41.6/40 

NVSS J202655+405408 

src3 

0 79 + 0 - 59 
u - ' y —0.45 

- 1 . 0 (> - 2 . 2 ) 

U.OO_o 36 

1 (fixed) 

0.19 X 10 -13 

1.27 X 10 “ 13 

28.1/25 

NVSS J202723+405706 

src4 

2.2(< 6.7) 

1 5 + 1 - 6 
10 -1.3 

- 

- 

0.06 X 10 -13 

1.47 X 10 “ 13 

16.0/25 

TYC 3156-1302-1 

src5 

q+15 

y — 7 

1.6(< 4.1) 

- 

- 

< 0.01 X 10 -13 

1.75 X 10 “ 13 

32.1/26 

NVSS J202642+405138 

src 6 

0(< 0.07) 

1 yo+O- 4 ^ 

A - ' -0.43 

- 

- 

0.48 X 10 -13 

0.86 X 10 “ 13 

18.1/13 

- 

Source 2 









srcl 

0 oo + 0 . 8 t> 

z.oo_o eg 

1 co+0.34 

1,0 -0.30 

- 

- 

0.24 x 10 “ 13 

5.46 X 10 “ 13 

120 . 8/101 

- 

src 2 

2.47(flxed b ) 

1 qo+0.73 
-0.83 

- 

- 

0.02 x 10“ 13 

0.92 x 10 “ 13 
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r is the photon index of the power-law model and kT and A are the temperature and abundance, respectively, in the 
plasma model. Abundance was fixed at the solar value if the value is not well constrained by a fit. / 0 . 5-2 and / 2-10 
are absorption-uncorrected fluxes in 0.5-2 keV and 2-10 keV, respectively. 


Objects named NVSS and TYC are radio sources ( Condon et al.l 1998ll and st ars I Hog et al. 1998l l. respectively. 


IRXS denotes a source detected in the ROSAT All-sky survey ( Voces et al. 1999l l. 


b N(}T) is not well determined and is fixed to the estimated Galactic absorption. 
c Crowded region with several possible counterparts (stars, radio sources and HII regions.) 
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4. Discussion 


4.1. Constraint on X-Ray Emission from the Cocoon 


As shown in Table 2 and Figure 4, no significant excess has been found in the X-ray spectra 
of the extended emission in any of the four regions. However, the data have been fitted with the 
CXB+GRXE(plus local diffuse X-ray emission) model described in § 3.2 and the fitted GRXE com¬ 
ponent may include part of the cocoon emission, if any. To search for an excess above the expected 
GRXE profile, we have examined the positional dependence of the CXB-subtracted intensity. We 
have focused on the hard-band (2-10 keV), in which the relative contribution of an extended X-ray 
emission from the cocoon is expected to be highest. 


Figure 5 shows how the intensity varies with Galactic latitude. The intensity of the CXB 
contribution is somewhat uncertain, due to uncertainty in the absorption and the field-to-field 
fluctuations of the CXB. To gauge the uncertainty of the absorption, we repeated the calculation 
of 77(H) with spin temperature of 100 K and optical-thin approximation, and found the impact on 
2-10 keV intensity is < 0.46x 10 -8 erg s _1 cm -2 sr -1 . According to lKushino et all (2002), the field- 
to-field fluctuation of the CXB in the AS'C'A-GIS field-of-view (~ 0.4 deg 2 ) is 5.1% (lex). Assuming 
that the fluctuation is Poisson’s noise, the expected fluctuation over the entire XIS field-of-view 
(~ 17.'8 x 17.'8) is 5.1 x a/0.4/0.088 = 11%. Therefore the uncertainty of the CXB contribution in 
the 2-10 keV range is (0.52-0.54) x 10~ 8 erg s^ 1 cm~ 2 sr -1 . U ncertainty in the NXB model also 
affects the intensity of the extended emission. iTawa et ajj (1200811 reported that the reproducibility 
of the XIS NXB model (lex) is ~ 5% and 3-4% in 1-7 keV and 5-12 keV ranges, respectively. 
We raised and lowered the NXB model by 5%, and found that the NXB-subtracted intensity in 
2-10 keV was altered by (0.26-0.63) x 10~ 8 erg s" 1 cm -2 sr^ 1 . We added the uncertainties due to 
the fluctuation of the CXB intensity and the NXB model reproducibility in quadrature, and took 
the effect from the assumption of spin temperature as systematic uncertainty. Even allowing for 
these uncertainties, the intensity profile (Figure 5) shows a clear monotonic decrease with increasing 
latitude, as expected for the GRXE. We have therefore concluded that most of the extended emission 
in the CXB-subtracted spectra in 2-10 keV is contributed by GRXE. 


We therefore impose constraints on the extended X-ray emission from the cocoon. Among 
the four positions, the Galactic latitude is the highest in Background 2. Therefore, we expect 
the GRXE intensity to be minimized there and take the intensity of Background 2 as the lower 
limit of the GRXE contribution in the source positions. A robust upper limit on the extended 
X-ray emission from the cocoon can be obtained by subtracting the intensity in Background 2. 
In this way, the upper limit in the 2-10 keV range for Source 1 and Source 2 was determined as 
3.4 x 10 -8 erg s -1 cm~ 2 sr" 1 and 1.3 x 10" 8 erg s -1 cm -2 sr -1 , respectively. 


We have also calculated the expected synchrotron X-ray intensity from a CR electron spectrum 
which can reproduce the GeV emission of the coco on bv up scattering the interstellar and stellar 
radiation fields modelled in the region (lAckermann et ah 201 lj) . The latter include the bright stars 
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of the Cyg OB2 and NGC 6910 associations, an average spectrum for the field stars, the intense 
IR field from the heated dust in the region, and the cosmological microwave background. The 
electron spectrum is based on the LIS one, which cuts off near 1 TeV, with power-law extensions 
up to 5, 10, 50, or 100 TeV. These spectra being too soft, we have applied an amplification factor 
of 120 x (E/50 GeV ) 0 ' * 5 to account for the Fermi -LAT data. It yields an electron density of 2.6 x 


10 " 9 m -3 TeV -1 at 1 TeV. We have assumed a thickness of 25 pc for the emission region along 
the line-of-sight (corresponding to 1° at d = 1.4 kpc in the plane of the sky), and a magnetic 
field strength of 20 /rG (deduced from pressure balance between the magnetic field and the gasj§. 
The synchrotron intensity expected in the 2-10 keV range from the electron spectrum extended to 
100 TeV is ~ 7.0 x 10 " 8 erg s " 1 cm -2 sr" 1 . As the upper limits from the Suzaku-XIS observations 
are 1/5-1/2 of this value, we infer that our data are inconsistent with a pure inverse-Compton 
scenario in 7 rays without a cut off in the electron spectrum. 

To further constrain the CR properties, we have investigated the multi-wavelength spectra, in¬ 
cluding radio and TeV 7 -rays as summarized in Figure El To compensate for the small field-of-view 
of the X-ray data compared to the size of the cocoon, we have integrated the averaged X-ray upper- 

, or 9.13xl0" 3 MeV 2 s " 1 cm " 2 sr " 1 MeV 


erg s 1 cm 2 sr 1 


limit (2-10 keV intensity of 2.35x 10 
assu min g a p hoto n index of 2.0) over the solid angle subtended by the cocoon (4.38 x 10 " 3 sr 


r -1 


Ackermann et ahl 12011). A potential TeV cou nterpart to the cocoon, MGRO J2031+41, has been 
found by MI LAGR O (Abdo et al. 2007 . 2012 '). More recent results from the ARGO-YBJ experi¬ 
ment (IBartoli et alJ 20141 ) confirm the presence of an extended TeV source of compatible position 
and size, even after subtracting the contribution from the known overlapping/ nearby TeV sources. 
Its spectrum smoothly connects to the Fermi -LAT data points as described by IBartoli et al.l (120141) 
(see also Fig. EIlJ- The radio flux has been averaged over the cocoon from the CGPS survey 
( Tavlor et al.l 12003 ). The values should be considered as upper limits to the possible synchrotron 
emission from the cocoon because free-free emission from the hot ionized gas of Cygnus X has not 
been removed. The WMAP free-free map is too coarse and Ha intensities are too absorbed in 
this region to allow a reliable subtraction of the free-free emission. The plotted curves for the syn¬ 
chrotron and inverse Compton models assume the amplified CR electron spectra described above, 
with the LIS cutting off near 1 TeV, and with power-law extensions to 5 TeV, 10 TeV, 50 TeV, and 
100 TeV cut-offs. Our X-ray data require a cut-off in electr on energy below 50 TeV, in agreemen t 
with the marked decline in TeV emission seen above 2 TeV (jAbdo et al.l 120121 : IBartoli et al.lI2014J ) . 
The X-ray data should therefore be regarded as an independent confirmation of the spectral cutoff 
in the electron scenario inferred from TeV observations. 


_ 6 7 Although this val ue is larger than the local interstellar magnetic field strength, we have adopted it as 

Ackermann et al] (2011) examined possible solutions to explain the GeV/TeV 7 -ray data while not violating the 

radio data in the specific environment of the cocoon. 

7 Although VERITAS has surveyed the cocoon region, it was not able to detect the 7 -ray emission on the angular 

sc ale of the Cygnus cocoon and MGRO J2031+41 with a standard ring-background estimation method, as described 


by Aliu et al. ( 20131 1 . 
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4.2. Properties of GRXE and Local Diffuse X-rays 


Finally we comment on the properties of GRXE and local diffuse X-ray emission. All four 
observations exhibit emission lines in the 0.5-0.6 keV range (see Figure 4). As shown in Table 2, 
the line central energy is ~ 0.55 keV and is not compatible with the O vn resonance line energy 
(0.574 keV) expected for the local hot plasma, but it is marginally consistent with the O vn forb id' 


den l i ne energy (0.561 keV ) which is a signature of a charge-exchange process (e.g.. lEbisawa et al 


20081 : F uiimoto e t al. 120071 b We therefore conclude that SWCX significantly contributes to our 
data below 1 keV. 

The intensity of the mid-temperature component (/ m id in Table 2) shows a positive correlation 
with Galactic latitude. This unusual behavi or for G RXE ma y be a ttributed to the structure of the 
Cygnus Super Bubble (CSB). As shown by U vaniker e t al. (120011 ). the positions of the Source 2 
and Background 2 images are inside or very close to the rim of the CSB, as seen by ROSAT in 
the 0.75 and 1.5 keV bands. Hence the intensity of the mid-temperature component should be 
significantly increased. We note that the -/V(H) m ;d and kT rn \ ( \ values of Source 2 and Background 2 
are not com patible with t he -/V (H) ~ 0.3 x 10 22 cm~ 2 and kT ~ 0.3 keV values ubiquitously found 
in CSB bv lKimura e t al. (|2013l). In particular, the Background 2 spectrum shows strong 0.7-1 keV 
emission (Fe-L and Ne lines), but lacks the line of Oviii, and therefore is not compatible with 
plas ma em is sion at ~ 0.3 keV. Although Suzaku provides high-quality spectra compared to those 
in Ki mura et al. (2013j) in terms of energy range and resolution, our data correspond to a much 


smaller region, so a detailed discussion of the origin of this discrepancy is difficult and beyond the 
scope of this study. 

The latitude dependence of the high-temperature component (/high in Table 2) gives a scale 
height of ~ l.°5, which is 2-3 times larger th an the GR XE scale height in the Scutum arm region 
close to the Galactic plane (~ 0.°5; e.g., Kaneda et al . 19971 1. This difference may be attributed 
to the proximity of Cygnus X (distance ~ 1.4 kpc), or to the dif ferent Ga l actic latitude ranges of 


the two studies (6 = 0.7-2.2 deg in our observations whereas iKaneda et al. ( 19971 ) analyzed data in 


b = 0-1.7 deg). The observed temperature (2-3 keV) is significantly lower than the value of hard- 


(Kaneda et al. 

1997 

) and l ~ 18° ( 

— i — rv i ■ ~ ^ ___ 

Uchivama et al. 

2009b) 


and ~ 5 keV were reported at l ~ 28° 
respectively. These Galactic longitude dependence of the temperature may provide useful informa¬ 
tion on the origin of GRXE (e.g., different populations of sources if GXRE is due to unresolved 
sources), but we defer a detailed discussion on this topic. 


5. Summary 

We have conducted a series of deep X-ray observations of the nearby star-forming region 
of Cygnus X using Suzaku XIS, in order to better understand the origin of the GeV emission 
revealed by Fermi -LAT in the cocoon of young cosmic rays pervading Cygnus X. After excluding 
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Fig. 5.— Dependence of the hard-band intensity (2-10 keV) on the Galactic latitude after subtract¬ 
ing the CXB contribution. The field-to-field fluctuation of the CXB intensity and the NXB model 
reproducibility are indicated by error bars, and the systematic uncertainty due to the absorption 
of the CXB model is shown by a shaded area. 

point sources and small-scale structures from the X-ray images, and after subtracting the CXB, 
we find that the X-ray distribution in 2-10 keV band monotonically decreases with increasing 
Galactic latitude, thus indicating that most of the extended emission relates to the GRXE. We 
have obtained robust upper limits to the diffuse X-ray emission from the 7 -ray cocoon. The results 
are incompatible with a pure inverse-Compton scenario for the origin of the 7 -ray emission without 
a cut-off in the electron spectrum above about 50 TeV, in agreement with the drop observed in 
the 7 -ray spectrum above 2 TeV. We have found a strong contribution of SWCX below 1 keV in 
our data. The mid-temperature component (kT = 0.6-0.7 keV) is likely to be significantly affected 
by the rim of CSB. The scale-height of the high-temperature component intensity (~ l.°5) is 2-3 
times larger than the scale height of GRXE in the Scutum arm region at l ~ 28°. We attribute 
this increased scale height to the proximity of Cygnus X. The observed temperature of the high- 
temperature component is 2-3 keV and is lower than found for the GRXE at smaller Galactic 
longitudes. 

We would like to thank K. Ebisawa, K. Hayashida and S. Yamauchi for valuable comments. 
We also thank the Suzaku team members for their dedicated support of the satellite operation and 
calibration. 
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Fig. 6.— Multi-wavelength data and models of the emission from the C ygnus cocoon regio n: the 


data points in the 0.1-100 GeV band come from the Fermi -LAT data f Ackerman n et alj 2011) 


t he d ata points and the p ower-law m ode l at TeY energies from the ARGO-YBJ (jBartoli et al 


20141 ) and the MILAGRO (jAbdo et ah 120121 ) experiments, respectively, the X-ray upp er li mit from 


t he p resent work, and the radio upper limits (triangles) from the CGPS survey (Taylor et_al 


20031 ). At TeV energies the MILAGRO data likely include all so urces aroun d the cocoo n area, but 


below 10 Te V the flux excee ds the sum from TeV J2032+4130 ( Aharonian et al.l 120051 ) and VER 
J2019+407 (Aliu et ah 20131 ) by at least a factor of two. Thin lines show the modelled synchrotron 
and inverse Compton emissions, assuming an amplified LIS electron spectrum, or this spectrum 
with power-law extensions to cut-off energies of 5, 10, 50, and 100 TeV. See text for details. 
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